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PACS 05 . 45 . Df - Fractals 

PACS 71 .30.+h - Metal- insulator transitions and other electronic transitions 

Abstract - Based on heuristic arguments we conjecture that an intimate relation exists between 
the eigenfunction multifractal dimensions D q of the eigenstates of critical random matrix ensembles 
D q i ~ qD g [q' + (q—q')D q ]~ 1 , 1 < q < 2. We verify this relation by extensive numerical calculations. 
We also demonstrate that the level compressibility \ describing level correlations can be related 
to D q in a unified way as D q = (1 — x)\l + (q — l)x] _ , thus generalizing existing relations with 
relevance to the disorder driven Anderson-transition. 
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Introduction. — It is well-known that the spatial 
fluctuations of the eigenstates in a disordered system at 
the Anderson-transition show multifractal characteristics 
|TJ[2] which has been demonstrated recently in a series of 
experiments [3]. Therefore the modeling and analysis of 
multifractal states has become of central importance pro- 
ducing many interesting results. For this purpose random 
matrix models have been invoked and studied recently 

(mi. 

Since the exact, analytical prediction of the multifractal 
dimensions of the states for the experimentally relevant 
Anderson-transition in d = 3 or the integer quantum-Hall 
transition in d = 2 seems to be out of reach, it is desirable 
to search for heuristic relations in order to understand the 
complexity of the states at criticality. In the present paper 
wc propose such heuristic relations that are numerically 
verified using various ensembles of random matrices. 

The spatial fluctuations of the eigenstates can be de- 
scribed by a set of multifractal dimensions D q defined by 
the scaling of the inverse mean eigenfunction participation 
numbers with the system size A^: 




always feel the entire system, i.e. D q — > d for all q. Multi- 
fractal states, on the other hand, should be described by 
the series of the D q , which are a nonlinear function of the 
parameter q. 

Spectral fluctuations can be characterized in many 
ways. A usual, often employed quantity is the level com- 
pressibility x, which is extracted from the limiting be- 
havior of the spectral number variance as T,^(E) = 
(n(E) 2 ) — (n(E)) ~ %E, where n(E) is the number of 
eigenstates in an interval of length E. The spectral fluc- 
tuations in a metallic system with extended states yield a 
vanishing compressibility, x ~* 0) while in a strongly dis- 
ordered insulating system the levels are uncorrelated, so 
they are easily compressible, \ = 1. However, for the mul- 
tifractal states an intermediate statistics exists, < x < 1' 
furthermore the spectral and eigenstate statistics are sup- 
posed to be coupled, which has been pointed out first in 
Ref. [7]. 

One of the most important generalized dimensions often 
used in this context is the information dimension D±. It 
is defined through the scaling of the mean eigenfunction 
entropy with the logarithm of the system size: 



N -(q-l)D q 



(1) 
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where (• • •) is the average over some eigenvalue window 
and over random realizations of the matrix. For strongly 
localized eigenstates these quantities do not scale with sys- 



tem size, i.e. D„ 



for all q, while extended states 



A further, well-known and widely used dimension is called 
the correlation dimension £>2, which is extracted from the 
inverse participation number from Eq. ([!]) using q = 2. 
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Fig. 1: D q as 
with 13 = 1. 
with Eq. (JSj) . 
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a function of b for the PBRM model at criticality 
The dashed lines are fits of the numerical data 



In a recent work [3] Bogomolny and Giraud have shown 
that in a d-dimensional critical system the information 
dimension D\ and the level compressibility \ are simply 
related as 

X + D 1 /d=l, (3) 

furthermore the generalized dimensions D q can be ex- 
pressed as 



D„ 



r(g-i/2) 
v^r( 5 ) 

1 ~QX > 



(i - x) , i - X « i 



(4) 



These expressions have been shown to be valid for various 
critical random matrix ensembles in Ref. [3]. 

As for the critical, three-dimensional Anderson transi- 
tion and the two-dimensional quantum-Hall transition it 
has been shown earlier that another relation holds between 
the level compressibility x an d the correlation dimension 
D 2 0: 

2x + D 2 /d = 1 . (5) 

This relation should obviously hold approximately only 
since < D^jd < 1 but < x < 1, leaving the range of 
validity for the limit of weak-multifractality. 

In the present work we show a series of relations be- 
tween various generalized dimensions, D q and D q i, and 
the level compressibility \ allowing for a generalization 
that for particular cases yields Eq. @j exactly and Eq. 
in the appropriate limit. In order to prove that, numerical 
simulations of various critical random matrix ensembles 
will be used. Further implications and more details will 
be presented elsewhere [5]. 

Model and heuristic relations. In Ref. 0] 
Eqs. <j3j) and (j4]) were shown to be correct numerically for 
the Power-Law Banded Random Matrix (PBRM) model 
2 . 9 , TU] at criticality. Below we will make use of this model 
to derive our main results. 

The PBRM model describes one-dimensional (Id) sam- 
ples of length N with random long-range hoppings. This 




Fig. 2: a q and j q = ai/a q as a function of q for the PBRM 
model at criticality with /3 = 1. The red dashed line equal to 
q is plotted to guide the eye. The error bars are the rms error 
of the fittings. 



model is represented by N x N real symmetric {j3 = 1) or 
complex hcrmitian (/3 = 2) matrices whose elements are 
statistically independent random variables drawn from a 
normal distribution with zero mean and a variance given 



by (\H„ 



(3- 1 and 



2 1 + [sin (n\m - n\/N) /{irb/N)] 



(6) 



where b and fj, are parameters. In Eq. ((5J) the PBRM 
model is in its periodic version; i.e. the Id sample is in 
a ring geometry. Theoretical considerations [21I5HH] and 
detailed numerical investigations [HH2JH2] have verified 
that the PBRM model undergoes a transition at fi = 1 
from localized states for fj, > 1 to delocalized states for 
fj, < 1. This transition shows all the key features of the 
disorder driven Anderson metal- insulator transition [2] , in- 
cluding multifractality of eigenfunctions and non-trivial 
spectral statistics. Thus the PBRM model possesses a 
line of critical points b G (0,oo) in the case of /i = 1. 
In the following we will focus on the PBRM model at 
criticality, fx = 1. By tuning the parameter b the states 
cross over from the nature of weak-multifractality {b ^> 1) 
which corresponds to extended-like or metallic-like states 
to strong-multifractality (d < 1) showing rather local- 
ized, i.e. insulator-like states. Meanwhile at the true, 
Anderson transition in d = 3 or at the integer quantum- 
Hall transition in d = 2, the states belong to the weakly 
multifractal regime, the PBRM model allows for an inves- 
tigation without such a limitation. The evolution of the 
generalized dimensions as a function of the parameter b 
therefore represent this behavior, i.e. D q — > 1 for i > 1 
and in the other limit of b 1 the multifractal dimensions 
vanish as D q ~ b PlfTU], 

Previously, for the PBRM model at criticality with 
P = 1, we have observed that both, D\ and D 2 can be 
approximated simply as Q3] D\ w [1 + (aife)^ 1 ] -1 and 
D 2 ~ [1 + (a^fr) -1 ] -1 where a\ y2 are fitting constants. 
This continuous function is a trivial interpolation between 
the limiting cases of low-6 and large-6 taking the half of 
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Fig. 4: (a-b) D q (open symbols) and Di[<? + (1 - g)Oi] _1 (full 
red symbols) [see Eq. (|13|l ] and (c-d) D q i (open symbols) and 
qD q [q + (q - q^Dq]^ 1 (full red symbols) [see Eq. l[I3)l] as a 
function of b for the PBRM model. 



3: (a) (1 - D q )[l + (7, - l)D q ] _1 and (b) (1 - D q )[l + analytically 0H] in the small- and large-6 limits: 
1)D 9 ] _1 [see Eqs. (TTT)) and (TT2)) ] as a function of b for 



Fig. 
(1- 

the PBRM model. The red dashed lines are the analytical 
prediction for x given in Eq. (flO|l . Inset in (b): qD q (l — Dq)^ 1 
as a function of 6, see Eq. (| lCjf) . The red dashed line equal to 
aib is plotted to guide the eye. 



the harmonic mean of the two as 
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Do 



= 1 



a q b ' 



(7) 



valid for q = 1 and 2. Here we generalize and propose 
the following heuristic expression for a wider range of the 
parameter q 



D„ 



[1 + (a q l 



(8) 



as a global fit for the multifractal dimensions D q of the 
PBRM model in both symmetries, (3 = 1 and j3 = 2. In 
Fig. [1] we show fits of Eq. ((8]) to numerically obtained D q 
as a function of b for some values of q and in Fig. [5] we 
plot the values of a q extracted from the fittings j^J We 
observe that Eq. ((5]) fits reasonably well the numerical 
D q for q > 1/2. It is important to stress that Eq. @ 
reproduces well the 6-dcpcndcncies predicted analytically 
[2] for the limits 6 <C 1 and b ^> 1. We noticed that by 
the use of Eq. ©, Eq. © leads to 



X « (1 + 



(9) 



which also reproduces well the 6-dependencies predicted 



The multifractal dimensions D ? were extracted from the linear 
fit of the logarithm of the inverse mean eigenfunction participation 
numbers versus the logarithm of iV, see Eq. JTJ. D\ was extracted 
from the linear fit of the mean eigenfunction entropy versus the log- 
arithm of N, see Eq. J2). We used N = 2", 8 < n < 13. The average 
was performed over 1 n ~ i eigenvectors with eigenvalues around the 
band center with 2 16— n realizations of the random matrices. 



1-46 
(27T&)- 1 



6< 1 
6> 1 



Then, by equating b in Eqs. (JSJ and © we get 
X w(l-D,)[l + (7,-l) J D,]- 1 , 



(10) 



(11) 



with 7 g = ttxjdq. We observed that 7 g w q in the range 
0.8 < q < 2.5, see Fig. [2J so in this range of q values we 
can write simplified relations between x and D q : 



I -D 



1 + {q - l)D q 



and D„ 



i-x 



1)X 



(12) 



The expression for D q in Eq. (fT2"]) reproduces Eq. (jl]) 

exactly for g = 1 and q = 2 and approximately for 

1 < g < 2.5. Moreover, Eq. (fUJ combined with Eq. © 
allows us to express any D g in terms of D\: 



D q ^D 1 [q+(l-q)D 1 ] 



(13) 



We also noticed that by equating x f° r different D q s 
form Eq. (fl~2"|) we could get recursive relations for them: 



9' 



1 - 1 - D, 



which lead to 13, 



9+1 



and D„/ = 



9^(1 



qD q 



q' + (q- q')D q 



(14) 



Dq) 1 , when 



g+1. These expressions also provide a relation between the 
correlation dimension and the information dimension or 
between the correlation dimension and the compressibility 
of the spectrum: 



Da = Di (2 - D l ) =(l-x){l + xY 



(15) 



It is relevant to add that in the weak multifractal regime, 
i.e. when x ~ * lj Eq. (fT5|) reproduces the relation given in 
Eq. © with d = 1, reported in [7]. 
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Fig. 5: Dx, D2, and Z?i(2 — Di) -1 as a function of a for the 
RSE (black open, red open, and blue full circles). Black and 
red dashed lines are the theoretical predictions for D± and D2, 
respectively, given in Eqs. (|20[) and (|21l) . The blue dashed line 
is the prediction for D2 given by Eq. (|22l) . 



Numerical results for the PBRM model. Here 
we verify the expressions (|1 Hll5f> for the PBRM model at 
criticality. Below we concentrate on the case (3 = 1 but 
we have already validated our results for /3 = 2. 

In Fig. [3] we plot (1 - D q )[l + (7, - 1)Z>J _1 and (1 - 
D q )[\ + (q — as a function of b for several values 

of q and observe good correspondence with the analytical 
prediction for \) that is, we verify the validity of Eqs. ([lip 
and (fT2"|) , respectively. In the inset of Fig. [SJb) we plot 



qD q (l — D q ) 1 as a function of 6, see Eq. ([T4"]) . which for 
the PBRM model acquires the simple form 

qD q (1 — Dq) -1 ps ga g 6 PS . 



(16) 



Then, in Fig. 2] we compare D q and Z? 9 ' with Z?i [q+ (1 — 
g)!?!]" 1 and qD q [q' +{q—q')D q ]~ 1 . respectively, for several 
values of q; that is, we verify the validity of Eqs. ([T3"]) and 
(fT4|) . Eq. (15]) is also validated in Fig.gjb). 

Additionally, in [15] the duality relation 



D 2 {B) + D 2 (B~ 



1 , £ = 2 1 / 4 tt6 , (17) 



was shown to be valid (with maximum deviations of 1%) 
for the PBRM model at criticality. We also want to com- 
ment that by the use of Eq. (JU) we could write D 2 (B) ps 
[1 + (SB)- 1 }- 1 , with 5 = a 2 /{2 1/ \), so relation (HI) gets 
the form 



D 2 (B) + D 2 (B~ 1 ) ps 1 — 



B(S-1) 



B + 5(B 2 + SB + 1) 



(18) 



We notice that the quantity D 2 (B) + D 2 (B~ 1 ) is very 
sensitive to the value of a 2 . So, the error in a 2 is magnified 
in the r.h.s. of Eq. (|TS|) . The maximal deviation from 1 
(of 7.3% and 2% for (3 = 1 and f3 = 2, respectively) occurs 
at B = 1 where the r.h.s. of Eq. (fT51) acquires the form 
1-[(S-1)/(S + 1)}\ 

Other critical ensembles. Remember that re- 
lations (|121I15[) were obtained form the combination of 
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Fig. 6: D q as a function of q for the RSE (black open symbols) 
for several values of a. Blue dashed lines are Eq. (|23|l . The 
green dashed line is Eq. (|24p with k — 2. The red dashed line 
is Eq. ([13]) with Zh(a = 1.5) = 0.96. 



Eqs. © and ©. That is, relations (|12j|15j) are expected 
to work in particular for the PBRM model at criticality. 
However, Eqs. (fT2"j) reproduce Eqs. and which were 
shown to be valid for the PBRM model but also for other 
critical ensembles [][] . Then the question is to which extent 
relations (|12M15[) are valid for critical ensembles different 
to the PBRM model. So, in the following we verify the 
validity of Eqs. (|12H15[) for other critical ensembles|j 

The Ruijsenaars-Schneider Ensemble (RSE). The 
RSE proposed in [TB] is defined as matrices of the form 



H mn = exp(i$ m ) 



1 — exp(27ria) 



iV[l-exp(27ri(m-n + a)/iV)] ' 



(19) 



where 1 < m < n, $ m are independent random phases 
distributed between and 27r, and a is a free parameter 
independent on N. When < a < 1, the compressibility 
and the multifractal dimensions take the form 0] 



\ 



(20) 



(a-iy and D q = 1 — q(a — 1) ; 
while in the vicinity of an integer k > 2, when \a — k\ <C 1, 
(a - kf/k 2 and D q = 1 - q(a - k) 2 /k 2 



X 



(21) 



As shown in g], Eqs. and (J5TJ) satisfy relation (0). 
Moreover, by direct substitution of Eqs. (j2T)|) [or Eqs. ([2T])] 
we verified that Eqs. (|12H15j) are also satisfied at leading 
order in (a — l) 2 [(a — k) 2 ]. 

In Fig. [5] we plot D\ and D 2 as a function of a for 
the RSE. Black and red dashed lines are the theoretical 
predictions for Di and D 2 , respectively, given in Eqs. (|2"U|) 
and (|2~T|) . As it was earlier shown in Ref. [4] , the analytical 
form of D q given in Eqs. ((20)) and ([21]) reproduces very well 



2 The multifractal dimensions i?i and D q for those ensembles 
were extracted numerically by the use of the same matrix sizes and 
ensemble realizations as for the PBRM model, if not indicated oth- 
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Fig. 7: Di, D2, and Di(2 — -Di) -1 as a function of g for the 
IQM model (black open, red open, and blue full circles). Black 
and red dashed lines are the theoretical predictions for Di and 
D2, respectively, given in Eqs. ()25|) . The blue dashed line is 
the prediction for D2 given by Eq. 



the numerically obtained D\. However, we notice that 
Eq. ([20]) does not describe well the numerical D2, mainly 
when a — > 0. Now, note that by plotting the numerically 
obtained D\/{2 — D±) we get good agreement with the 
numerical data for D2 , that is Eq. (fT5)) works well for this 
model. Then, if we take D\ w 1 — (a — l) 2 and D\ w 
1 — (a — k) 2 /k 2 as theoretical predictions for D\ and plug 
them into Eq. (|T5|) we get 



Do 



l-(«-l) 2 
l + (a-l) 2 



and D2 



kf 



k 2 + (a-ky 



(22) 



for < a < 1 and |o — k\ <C 1 with k > 2, respectively; 
which in fact work much better than D2 ~ 1 — 2(a — l) 2 
and D2 ~ 1 — 2(a — k) 2 /k 2 , correspondingly; see Fig. [5J 

To get expressions for D q we substituted x ~ (a — I) 2 
and x ~ (a - A:) 2 / fc2 [ or f 1 w 1 - (o - l) 2 and L»i w 
1 - (a - fc) 2 /fc 2 ] into Eq. {HJ [or Eq. ^3])]. to get 



D q ^ [l-(a-l) 2 ] [l + (g-l)(a-l) 2 ]" 



and 



D q P3 [fc 2 - (a - fc) 2 ] [A: 2 + (g - l)(a - fc) 2 



(23) 



(24) 



In Fig. [6] we plot D q as a function of q for the RSE for 
several values of a. We also plot Eqs. (f2"3"]) and and 
observe rather good correspondence with the numerical 
data mainly in the range 1 < q < 2. Notice that neither 
Eq. d23I) nor Eq. (J53) can be used for a = 1.5. For that 
case we substituted the numerically obtained value of D\ 
into Eq. (|13p and again observe good correspondence for 
< q < 2, see the red dashed line in Fig. 

Intermediate quantum maps. A variant of the RSE 
was studied in |17j with the name of intermediate quantum 
maps (IQM) model. In this model the parameter a of 
the RSE equals cN/g with cN = ±1 mod g, being g the 
parameter of the IQM model. For the IQM model the 




Fig. 8: D q as a function of q for the IQM model (black open 
symbols) for g = 2, 3, 10, and 18. Blue dashed lines are 
Eq. (|27p . Red dashed lines are Eq. I|13|) where the numeri- 
cal values Di = 0.538, 0.683, 0.904, and 0.947 have been used 
for g = 2, 3, 10, and 18, respectively. Black dashed lines are 
D q from Eqs. 
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Fig. 9: Dq as a function of q for the CUE (taken from \18\). 
The blue dashed line is D q from Eq. ((29j) for J/W = 0.1. Red 
dashed lines are the prediction for D q given by Eq. (| 13p using 
Di = 0.2805 and 0.8443 for J/W = 0.1 and 3, respectively. 



compressibility and the multifractal dimensions take the 
form [IT] 



X 1/g and D q K,l- q/g . 



(25) 



As for the RSE, here Eqs. (J2SJ satisfy relation (0). 
Again, by direct substitution of Eqs. (|25p we verified that 
Eqs. (|124I15|) are satisfied at leading order in 1/g, j> 1. 

We want to mention that in [T7] it was shown that 
Eq. (f2"5"j) reproduces well the numerically obtained D\ but 
underestimates the numerical D2, in particular for small 
<7, see Fig. [71 Now, notice that by plotting the numerically 
obtained D\/(2 — D{) we nicely reproduce the numerical 
data for D2 , that is Eq. (fT5")) works well also for this model. 
Then, if we take D\ ~ 1 — 1/g as the theoretical prediction 
for Di and plug it into Eq. (jl~5j) we get 



Do 



(1 - 1/g) (1 + l/ g y 



which in fact works much better than D2 
reproducing the numerical D2, see Fig. [JJ 



(26) 
1 - 2/g in 
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To get the expression for D q we substituted \ ~ Vff or 
Di ps 1 — 1/g into Eq. (fT2"|) or (fT3|) . respectively, to get 



(<?-!)(<? + ?- 1)" 



(27) 



In Fig. [5] we plot D q as a function of q for the IQM model 
for some values of g. We also plot Eq. (|27p and observe 
that it falls below the numerical data mainly for small g. 
However, by substituting the numerically obtained values 
of Di into Eq. (fT3|) we get much better correspondence 
with the numerical D q , mainly for 1 < q < 2. In Fig. [5] 
we also include D q from Eq. (|2"5|) . We may conclude that 
while Eq. (|25p reproduces well the numerical for q < 1, 
Eq. (|2"7]l can serve as the analytical continuation for q > 1. 

The critical ultrametric ensemble. The critical ul- 
trametric ensemble (CUE) proposed in [T5] consists of 
2 K x 2 K Hermitian matrices whose matrix elements are 
Gaussian random variables with zero mean and variance 



(\H mm \ 2 )=W 2 , (\H mn \ 2 ) 



71 J 



(28) 



where d mn is the ultrametric distance between m and n 
on the binary tree with K levels and the root of 1. The 
parameter in this model is the ratio J/W. For the CUE, 
when J/W <§C 1, the compressibility and the multifractal 
dimensions have the form {4"1IT8"] 



X 



J 



and D n 



J 0Fr(g - 1/2) 
W v / 21n2r( g ) 



(29) 



Eqs. (2ni) satisfy relation © at first order in J/W g]. 
Again, as for the previous critical ensembles, by direct sub- 
stitution of Eqs. ((29)) we verified that Eqs. (|12II15[) are sat- 
isfied at leading order in J/W, for 0.8 < q < 2.5; because 
in this range of q we have that T(q — 0.5) /^T(q) ps 1/q, 
In Fig. IH1 we show D q as a function of q for the CUE 
for J/W = 0.1 and 3. The data was taken from [Tg] . 
The blue dashed line is D q from Eq. (g9) for J/W = 0.1. 
Notice that since Eq. (f2T)]) is only valid when J/W -C 1 
and for q > 3/4 one can not use it to predict D q for 
J/W = 3. However, with Eq. (fLl]) using as input the 
numerically obtained D\ we got good predictions for D q 
for small and large values of J/W and even for values of 
q smaller than 3/4. This is shown in Fig. [S] where we plot 
Eq. ([13]) (red dashed lines) using D 1 = 0.2805 and 0.8443 
for J/W = 0.1 and 3, respectively The values of D\ were 
obtained by the interpolation of the D q data. We observe 
good correspondence between Eq. (fl~3|) and the numerical 
D q for < q < 10. 

Conclusions. — In this paper we propose heuristic re- 
lations on one hand between the generalized multifractal 
dimensions, D q and D q i , for a relatively wide range of the 
parameters q and q' , and on the other hand between these 
dimensions and the level compressibility \- As a result we 
find a general framework embracing an earlier [7] and a 
recent one [J. Our proposed relations have been backed 
by numerical simulation on various random matrix ensem- 
bles whose eigenstates have multifractal properties. These 



results call for further theoretical as well as numerical in- 
vestigations. 
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